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Correlation functions of the open XXZ chain II 
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Abstract 

We derive compact multiple integral formulas for several physical spin cor- 
relation functions in the semi-infinite XXZ chain with a longitudinal boundary 
magnetic field. Our formulas follow from several effective re-summations of the 
multiple integral representation for the elementary blocks obtained in our pre- 
vious article (I). In the free fermion point we compute the local magnetization 
as well as the density of energy profiles. These quantities, in addition to their 
bulk behavior, exhibit Friedel type oscillations induced by the boundary; their 
amplitudes depend on the boundary magnetic field and decay algebraically in 
terms of the distance to the boundary. 

1 Introduction 

The Hamiltonian of the Heisenberg XXZ spin-1/2 finite chain [1] with diagonal 
boundary conditions (namely with longitudinal boundary magnetic fields) is defined 
as [2,3] 

M-l 

H =Y1 + + A ~ 1) } + h-af + h + a z M . (1.1) 

m=l 
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This is a linear operator acting in the quantum space !K = <8) 3^ m; — C 2 , of 

m=l 

dimension 2 M of the chain. In this expression, a^, &m denote local spin operators 
(acting as Pauli matrices) at site m, A is the anisotropy parameter and h± are the 
boundary (longitudinal) magnetic fields. 

We have recently developed a method to compute the so-called elementary blocks 
of correlation functions for this model (see [4], that we refer to as Paper I in the 
following) in the framework of the (algebraic) Bethe ansatz [5-18] for boundary in- 
tegrable systems [2,3,19-31]. The results essentially agree with previous expressions 
derived from the vertex operator approach [32,33]. The purpose of the present paper 
is to obtain the physical spin correlation functions for this model, in particular, the 
one point functions for the local spin operators at distance m from the boundary as 
well as several two point functions (like boundary-bulk correlation functions). There 
are numerous physical interests in such quantities that can be measured in actual 
experiments [34-46]. 

In much the same way as in the bulk case [47-52] , the computation of the physical 
correlation functions amounts to obtain effective re-summations of the multiple inte- 
gral representations derived for the elementary blocks. For example, the one point 
functions at distance m from the boundary, such as the local magnetization (c^), 
can be written as a sum of 2 m elementary blocks. We will show how to obtain com- 
pact expressions for such objects, typically involving the sum of only m terms, each 
containing multiple integrals whose integrants have a structure similar to the one of 
the elementary blocks. In the free fermion point we are able to compute these multi- 
ple integrals (and hence the corresponding correlations functions) almost completely 
by reducing them to single integrals. For instance the local magnetization and the 
density of energy profiles (a quantity of interest in the study and the understanding 
of the interplay between quantum entanglement and quantum criticality [53-61]) are 
expressed as single integrals. Hence, their asymptotic behavior at long distance m 
from the boundary can be explicitly evaluated. In addition to the bulk constant value 
they exhibit Friedel type oscillations [44-46,60,61], algebraically decaying with the 
distance m, their amplitudes being rational functions of the boundary magnetic field, 
in agreement with field theory predictions [38-42,62-70]. 

We start this paper with a short technical introduction concerning the algebraic 
Bethe ansatz approach to the open XXZ spin-1/2 chain subject to diagonal boundary 
magnetic fields. This preliminary section is followed in Section by a reminder 
of the method proposed in [4] to compute correlation functions of open integrable 
models in the framework of algebraic Bethe ansatz. In Section [5] we obtain formulae 
for the action of local operators on arbitrary boundary states in a form suitable for 
taking later on the thermodynamic limit. Using these results, we derive a series 
representation for the generating function (Q m (k)) of bulk-boundary a z correlation 
functions in Section This formula is the boundary analogue of the original series [49] 
in the bulk case. In Section^ we obtain a formula for {Q m (k)) alternative to the one 
inferred in Section [5j We also give multiple integral representations for an d 
for the local density of energy. These formulae are obtained by a direct resummation 
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of the corresponding elementary blocks. It is worth stressing that we actually have two 
representations for their integrand. The first one is in the spirit of the bulk case [52] 
and involves the Izergin determinant representation [71] for the partition function of 
the six vertex model with domain wall boundary conditions. The second one involves 
the Tsuchiya [72] determinant representation for the partition function of the six 
vertex model with reflecting ends. The next section is devoted to the free fermion 
point. For that case are able to reduce the multiple integrals to one dimensional ones. 
This allows us to write the leading asymptotics of the local magnetization and of 
the density of energy profiles as well as of the {&m+i (T i) correlation function. Our 
conclusions are presented in the last section. 



2 The open XXZ spin-1/2 chain 



The spectrum of H can be obtained by algebraic Bethe ansatz (ABA) [3]. The 
central tool of this method is the boundary monodromy matrix, which will be defined 
after we introduce some necessary notations. 

Here and in the following we adopt the standard parameterizations A = cosh 77 
and h± = sinh77coth£±. 

Let R : C — > End(V <g>V), V ~ C 2 , be the i?-matrix of the six- vertex model, 
obtained as the trigonometric solution of the Yang-Baxter equation: 



R(u) = sinh (u + 77) R(u) , with R(u) = 



/I 









b(u) 
c(u) 






c(u) 
b{u) 







1/ 



and 



b(u) 



sinhu 



sinh(n + 77) ' 
The bulk monodromy matrix T(A) 6 End(Vo 



c{u) 



sinhr/ 



(2.1) 



(2.2) 



sinh(ii + rj) 
"K), Vq ~ C 2 , is defined as an ordered 



product of R matrices: 

T (A) = Rqm (A — £m) ■ • ■ Rqi (A - 6) = 



A(X) B (A) 
C(A) D (A) 



(2.3) 



The subscript labels here the two-dimensional auxiliary space Vq, whereas subscripts 
m running from 1 to M refer to the quantum spaces Jf m of the chain. Besides, we 
attach an inhomogeneity parameter £ m to each site m of the chain. We recall that 
T (A) satisfies the Yang-Baxter algebra, on Vq (g) Vq> "K: 



Roq, (A - n) T (A) T G , Qu) = T (A) T v (ji) i? 00 ' (A - y) 



(2.4) 



Let us also introduce the two boundary matrices, K± (A) = K (A ± rj/2; £±), where 
K (A; is the 2x2 matrix acting on the auxiliary space: 



K(\;0 



sinh {X + 

sinh (f - A) 



[0] 



(2.5) 
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The boundary monodromy matrix U (A) [3] is built out of a product of T (A) and 
if ^ (A), namel}0 

W t0 = T tn (A) (A) ?o° (A) = ( £(A) p (A) ) ^ ' (2 - 6) 

where 

f (A) = 22 10 (A + Ci -??).. • #mo (A + £ M ~ rj) 

M 

= ("1) M II[ sinh ( A + 0) smh (A + - V)} T^ 1 (-A + r?) . (2.7) 

3=1 

This boundary monodromy matrix satisfies the reflection algebra first introduced 
in [20]: 

#00' ("A + n) U$ (A) R m , (—A — /i — r?) tifr (/*) 

= W V (/x) #oo' (-A - /x - rj) Ul° (A) iioo' (-A + n) . (2.8) 

The commuting charges of the XXZ spin- 1/2 chain with diagonal boundary con- 
ditions are realized by the one-parameter family of transfer matrices: 

T (A) = tr [ Uq (A) K- (A)] , (2.9) 

and the Hamiltonian (jl.ip is obtained in terms of the derivative — ^ - \ \ =r) /2 in the 

homogeneous case £j = rj/2, i = 1, . . . , M. 

Common eigenstates of all transfer matrices (and thus of the Hamiltonian in 
the homogeneous case) can be constructed by successive actions of B (A) operators on 
the reference state 10 ) which is the ferromagnetic state with all the spins up. More 
precisely, the state □ 

|{A}?) 6 = £(A 1 )...£(A n )|0) (2.10) 

is a common eigenstate of the transfer matrices if the set of spectral parameters 
{A}" = {Aj}i<j< n is a solution of the Bethe equations 

Vj (A i5 {A}?) = Vi (-A j ; {A}™) , j = 1, . . . , n, (2.11) 

where 

y j (x;{Xr i )= ^ ;{A} " ) , , 
s(Aj,x — 77) 

y (x; {A}™) = —a (x) d (— x) sinh (x + £ + — r?/2) sinh (x + £_ — 77/2) 

n 

X JJ«(x-»7,Ai) . (2.12) 



1 Note that it corresponds to the matrix U+ of our previous article (I). Since we consider only the 
'+' case in the present article, we do not specify it in the notations. 

2 In order to lighten the formulae we have slightly changed the notation with respect to the one in 
Paper I. Namely, the vector | {A}™ )& corresponds to \tp+ ({A}) ) in [4]. Such a boundary state should 
in particular be distinguished from the corresponding bulk state that we merely denote | {A}™ ). 
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Here and in the following, s (A, (i) denotes the function 

s (A, n) = sinh (A + /i) sinh (A - /x) , (2.13) 

and the functions a (A) and d (A) stand respectively for the eigenvalues of the bulk 
operators A(X) and -D(A) on the pseudo-vacuum |0): 

M M 

a(A) = fJsinh(A-ei+r/), d (A) = JJ sinh (A - &) . (2.14) 

i=l i=l 

Of course it is also possible to implement the Bethe ansatz starting from the dual 
state (0 | and acting on it with C (A) operators: 

6 ({A}? | = (0|C(Ai)...C(A n ). (2.15) 

The description of the ground state of Ti in the half-infinite chain depends on the 
regime. One should distinguish the two domains — 1 < A < 1 (massless regime) and 
A > 1 (massive regime): 

7T ~ 7T 

otj = Xj, ( = irj>0, £_ = with - — < £_ < — , for - 1 < A < 1, 

atj = C = -i) > 0, £- = -I- + »<*77, w ith £_ G R, for A > 1, 

where 5=1 for < sinh£ and <5 = otherwise. Thus, to a given set of roots {Xj} 
corresponds a set of variables {ay} given by the previous change of variables. Note 
that the nature of the ground state rapidities depends on the value of the boundary 
field h-. 

Indeed, when ^_ < or > £/2, the ground state of the Hamiltonian (|1.1|) is 
given in both regimes by the maximum number N of roots Xj corresponding to real 
(positive) a.j such that cosp(Aj) < A. In the thermodynamic limit M — > oo, these 
roots Xj form a dense distribution on an interval [0, A] of the real or imaginary axis. 
Their density 

p(Xj)= lim [MiXj+x-Xj)}- 1 (2.16) 

M— *oo 

satisfies the integral equation 

A 

f zsinh(2r?) , 2isinhr? 

2irp(X)+ / r~ — ^-p A dA = — (2.17) 

^ V 1 J S (A- M ,77)^ 1 ; s(A,r//2)' V ; 

-A 

with A = +oo in the massless regime, and A = —itt/2 in the massive one. The density 
can be expressed in terms of usual functions: 



P(A) 



' -1 < A < 1; 



C cosh (.A/C) ' , 

i „ /sinhnCA 3 (i\;-Q ( " " , 



^_ — ;^_vj^»_v/ 1 < A 
, vr „>! VcoshnC/ 9 4 (iX; -Q ' 
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However, when < < £/2, the ground state also admits a root A (corresponding 
to a complex a) which tends to rj/2 — with exponentially small corrections in the 
large M limit. In that case, the density of real roots is still given by the solution of 

3 The ABA approach to correlation functions 

A zero temperature correlation function is the normalized expectation value, in 
the ground state of the Hamiltonian (|l.ip . of some /oca|l operator O m , 



b( {A}f |O m |{A} f ) b 
><{A}f |{A}f 



(°™) = ^I^XnI (3-1) 



where the parameters A are the solutions of the ground state Bethe equations. 

In order to compute such a correlation function, one should first derive the ac- 
tion of the corresponding local operator on the boundary state | {A}^)?,, and then 
evaluate the resulting scalar products. We have constructed in [4] a method to solve 
this problem. This method is based on a revisited version of the quantum inverse 
problem, first introduced in [47, 73] for the XXZ spin chain with periodic boundary 
conditions. Once a local operator is reconstructed in terms of the entries of the bulk 
monodromy matrix, its action on boundary states can then be computed thanks to 
the decomposition of boundary states in terms of bulk states and to the Yang-Baxter 
commutation relations. 

We shall now recall the main points of our method. 

3.1 The bulk inverse problem revisited 

Proposition 3.1 (Solution of the bulk inverse problem) [47,73] Let E m be an 
elementary matrix acting non-trivially only on the m th site of the chain, then 

m—X m 

E m = II ( A + D ) (&) tr ( T o (U) E$) H(A + D)- 1 (6) . (3.2) 

k=l k=l 

Note that, thanks to the crossing symmetry of the i?-matrix, one can recast the 
inverse of the bulk transfer matrix at inhomogeneity parameter {A + D)^ 1 (£&) in 
terms of the transfer matrix at shifted parameter (A + D) (£& — rf), namely 

(A + g)- (fa) + f^ . (3.3) 
a (£fc) d - rj) 

It is worth pointing out that the products of elementary matrices on the first m 
sites of the chain define a basis in the space of local operators O m , so that (13. 2D allows 
one to define a reconstruction for all such operators. However, this reconstruction is 



„ in 

i.e. acting non-trivially only in (g) !Hfc. 
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especially convenient when one wants to obtain the action on a bulk Bethe state; 
indeed, in such a case, the product of bulk transfer matrices merely produces a nu- 
merical factor. This is no longer the case when one acts on a boundary Bethe state. 
The theorem below allows one to reconstruct local operators in a way adapted to an 
action on boundary states. 

Theorem 3.1 [4] For any set of inhomogeneity parameters ■ ■ ■ the prod- 

uct of bulk operators 

T «n < n &») • • • T tii , n (&) T- Hx % (& -ri)...T- Hn - <n - „) (3.4) 

vanishes if, for some k € . . . , i n }, = e*.. 

Thus we have : 

Corollary 3.1 A product of elementary matrices acting on the first m sites of the 
chain can be expressed as a single monomial in the entries of the bulk monodromy 
matrix: 

m 

i=i 

x T e , ei (6 ) . . . T eL em (U) T- em Sn (U-V)--- T- ei - ei (6 - V) (3-5) 

with Ei = e' i + 1 (mod 2). 

This result represents a strong simplification. Indeed, it means that, over the 
2 m monomials appearing in the reconstruction of a local operator (|3.2p . only one is 
non-zero. We shall now explain how to compute the action of this non-vanishing 
monomial on an arbitrary (bulk or boundary) state. 

3.2 Action on bulk and boundary states 

Before stating the lemma which explains how to derive the action of the former 
monomial on a bulk state, we recall that the action of A(n) or D(fi) on a bulk state 
IjA}^) = rii=i produces two kinds of terms: the direct term, where all 
rapidities remain unchanged, and indirect terms where one Xj is replaced by [i. 

Lemma 1 (Action on a bulk state) [4] The action on a bulk state \ {X}^ ) of a 
string of operators 

S v ' S v ' 

(1) (2) 

with Ei = e[ + 1 (mod 2), satisfies the restrictions: 

• The only non-zero contributions of the tail operators (2) come from 
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(i) the indirect action of all A(£i — rf) operators; 

(ii) the direct action of all D(£i — rf) operators. 

• In what concerns the head operators {!), 

(in) if e[ = 1, the action of the operator T e / ej (£i) (i.e. A(£i) or B(£i)) does not 
result in any substitution of a parameter £j — rj; 

(iv) ife'i = 2, the action of the operator T e i £i (£z) (i.e. D(£i) orC(£i)) substitutes 
£l — rj with fa; moreover, if there were others parameters £j — rj, j 7^ I, in 
the initial state, they are still present in the resulting state. 

This lemma enables us to compute the action of local operators on any (arbitrary) 
bulk state. In order to compute the action on a boundary state, we use the fact that 
the latter can be decomposed in terms of bulk states: 

Proposition 3.2 (Boundary-bulk decomposition) [4], [74] Let | {A}™ }& be an 
arbitrary boundary state, then it can be expressed in terms of bulk states as 

|{A}™>6= £ #f CTi} ({A}?)|{A CT }?>, (3.7) 



<Tj=± 

i=l,...,n 



with 



< } (W)=n<(A,) • n smh( 5f7 ) - ^ 

1 1^ J> smh(A r J 

j=\ l<r<s<n y rs ' 

In this expression, H®(\) denotes the "one-particle" boundary-bulk coefficient, which 
can be written as 

H*(\) = a (-l) M d(-X°) Sln ^ 2 h A 2 + V) sinh (A CT + - r//2) . (3.9) 
Here we have used the notations: 

Ars = X r — X s , X rs = X r + A s , (3.10) 



xy 



o~jXj, and more generally {A CT }™ = {aJ j }™ . (3.11) 



It is remarkable that, by using this decomposition and the previous lemma, we 
are able to express the action of a local operator O m on an arbitrary boundary state 
as a linear combination of such boundary states: 

O m \ {A}f } b = Y J <?a m ({A} ; {£}) I M ieam ) b , (3.12) 

Q-m 

where the summation is taken over certain subsets {l^i} iea of {A}^ U and 
where C am are coefficients which can be computed genericalbyE]. 



See section 5.3 of [4] for the explicit expression in the case of a product of elementary matrices. 



s 



3.3 Prom scalar products to correlation functions 

It now remains, in order to obtain the correlation function (13.11). to take the scalar 



product of this resulting combination of states with the ground state ^({A}^ |. This 
can be done by using the trigonometric generalization [4] of the rational [74] formula 
for the scalar product between a boundary Bethe state and an arbitrary boundary 
state. In particular, we have to evaluate the following type of renormalized scalar 
product: 

S({A},M) = b ^j|j^ )fc , (3.13) 

where the sets {A} and {//} are partitioned according to: 

{A} = {A a } aea _U{A 6 } 6eQ+ , {M} = {A«} oea _U{6} 667+ , (3.14) 

with \a + \ = |7+|. Here the parameters A are the solutions of the ground state bound- 
ary Bethe equations, {Cfe}fog 7+ are arbitrary inhomogeneities and a+Ua_ is a partition 
of {1,...,N}. 

Since we are especially interested in the thermodynamic limit M — > +oo of the 
correlation function (|3.1|) . we only recall the explicit formula for the leading asymp- 
totic contribution of the renormalized scalar product [4] : 

S ({ A} , M) = (Jj) ' + ' §({A W , {£} 7+ ; { A} Q _ ) det + [Sat] • (3.15) 

be-y+ 

The coefficient §({A} a+ , {^} 7+ ; {A} Q _ ) has been computed in [4]. Note that, since 
{A} is a solution of the boundary Bethe equations, there is some sign arbitrariness in 
the expression of this coefficient: indeed, it is in fact equal to 

11 s(Ah,A a ) rr \ \ 

S CT ({A} a+ ,{£} 7+ ;{A} Q _) = — — — 11 n 771 \ \ 

a>b 

n y(&, {A} a+Ua _ ) sinh(2&, + rj) yr sinh(2A^ - rj) sinh(2A a ) 
sinh(26) a |l + y(A-{A} Q+Ua _)sinh(2A a +r ? )' [ ' ) 

for any value of a a £ {+,—}, a 6 a + , where y is the function defined in (|2.12p . 
When M is large, the matrix elements of S reduce to 

~ f 2*7rMsinh(A a -£_ + »7/2)tf(A ,&) if A a = A , 

O r ,h ~ \ ~ (6. It) 

m-+oo \ p- 1 (A a ) * (A a , 6) if A a + A , 

the corrections being of order O (1/M), and 

'w '^Lg-T* - <3 - 18) 
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The determinant structure of the scalar product as well as peculiarities of the 
coefficients C am enable us to write: 

(°™) = J^ E ^(Mi,{e}D(l + 0(l/M)), (3.19) 

|/|=m 

in which the coefficient H m ({u}j , {£}™) can be computed generically. Taking the 
thermodynamic limit M — > +00 we recast the sums over replaced rapidities A into 
integrals: 

1 w /■ 

m E E * / ( A ^) M — ^ / dA p (A) ; (A) ' v/ G c ° m ■ (3 - 20) 

i=l <7i=± ^ °° i 

In the boundary model the contour of integration depends on the anisotropy pa- 
rameter A and on the boundary field 



4 Action of local operators on boundary states 

Using Corollary 13. 1|, Lemma Q] and the boundary-bulk decomposition of Propo- 
sition 13.21 it is easy to compute the action of a product of elementary matrices of 
the form (13. 5p on an arbitrary boundary state. This computation was explicitely 
performed in [4], and enabled us there to obtain some expressions for the elementary 
building blocks of correlation functions. 

The aim of the present article is to obtain such expressions for physical correlation 
functions, and in particular for one-point functions. Therefore, if we want to use the 
method recalled in Section [31 the main problem is to obtain some resummed formulas 
directly for the action of the local spin operators we consider. This is the purpose of 
the present section. 

In the first part of this section, we derive the action of the operator 

m mm 

Qm («) = [J {E} 1 + KEf ) = (A + kD) (£0 H(A + D)- 1 &) (4.1) 

i=l i=l i=l 

on arbitrary boundary states. {Q m (k)) can be interpreted in the boundary model as 
the generating function of the magnetization at a distance m from the boundary: 

( 1_^ L) = DMQm (k))U=1) ( 4. 2) 

where D m is the lattice derivative : D m u m = u m+ i — u m . 

Then, in the second part of this section, we give the formulas for the action of the 
local spin operators E^ = 1 ~^ m , E^ = cr+ and E^ = cr~ on arbitrary boundary 
states. Note that the action of E}} follows from the fact that E}} = 1 — E?3. 
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4.1 Action of Q m (k) 

We start by computing the action of Q rn (n) on an arbitrary bulk state, and then 
infer from this formula its action on arbitrary boundary states. 

Proposition 4.1 The action of Q m (k) on an arbitrary bulk state | {A}^ ) can be 

expressed as 

m 

Qm (*) i {A}f ) = E R n to . n) i {ei 7+ u w a _ > . (4.3) 

«=0 Pa ; ^ 



In i/te above formula, we sum over all possible partitions V\ and of the sets {A} 1 
and {^} m into subsets {A} a+ U {A} a and {£} 7+ U {£} 7 _ respectively, satisfying the 
constraint on the cardinality \ a + | = | 7+ |= n: 

Vx: {A}f = {A} Q+ U{A} a _, \a + \=n, (4.4) 

V(, : {£}? = tf} 7+ U «} 7 _ , |7+I=n- (4-5) 

T/ie coefficient (V\ , "Pg) spZife into too parts, 

KCPx,V i )=R(Vx,P^) ^(U} 7+ ,{A} Q+ ), (4.6) 

t/ie /jrst one having a product structure, 



N 



n «(aj n / (a&, A a ) n 
r (vx , p 5 ) = — ^ L n b u f(x rv (47) 

n Ufa) n <* 7 - 

£167+ L ftea_UQ! + J 

and the second one, which depends here only on the subsets {A} Q+ and {£} 7+ , being 
given as a ratio of two determinants, 

^(M?,M?) = det B [M < ,(MJ,MJ)] det- 1 ^,^^^^^ ) . (4.8) 

The entries of the matrix M K read 

f(Pa,fij) T T f(H^a) 

and the functions f and t stand for 



w(.^w;)]^««)-«'(«)n|Sn|s • (4 



.9) 



± , . \ sinh 77 , . . , sinh (A — /x + 77) , 1 

t(A ^ ) = sinh ( A-^)sinh(A-^ + n) ' 7 ( ^ } = sinh(A-/ - (410) 
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The above theorem appears as a non-trivial generalization of the action of Q m (k) 
on bulk Bethe eigenvectors [49]. Indeed, when | {A} ) is not an eigenstate of the bulk 
transfer matrix, then H2=i + D) (£j) does not act by multiplication any more. Of 
course our result reproduces the previous case when we send the parameters A to a 
solution of the bulk Bethe equations. 

Proof — The proof goes by induction on m. 

Property (|4.3p is obvious for m = 1. Assume that it holds for some m. To prove 
its validity for m + 1 we have to compute 

Q m+1 (k) \ {X\ } = — — \{\\ ). (4.11) 

a(£ m+ i)d(£ m+ i -77) 

Let us first reproduce the coefficient R 1 ^ (V\ t'Pf) i n t ne case when the partition 
V% is such that £ m +i {£} 7+ - The corresponding state | {A} Q U {£} 7+ ) can on ly De 
obtained by the direct action of (A + kD) (£ m +i). In order to reproduce the claimed 
form of the coefficient it is enough to prove that (A + kD) (£ m +i — 77) acts directly. 
Suppose that this is not the case. Then Q m (k) acts on a state containing £, m +i — f]- 
In virtue of Lemma [JJ the action of Q m (k) on these states cannot replace £m+i — V- 
Thus (A + D) (£ m +i) exchanges £ m +i — 77 with £ m +i, which leads to a contradiction. 

We still have to reproduce the coefficient R* (T\ , V^) corresponding to states 
I {-^}q_ U i?J"7+ ) sucn that Cm+i £ {?)7 + - Theorem 13.11 yields the decomposition: 

m+l) ^Cm V^) \ 1 f \ 

Q m+1 (k) = — — — D (£ m+1 - r?) +—7z rjTZ r A (£ m+ i - rj) , 

(1) (2) 

whereas Lemma [T] ensures that 

• (1) only acts directly; indeed A (£ m +i) cannot replace £ m +i — V by £ m +i; 

• (2) acts indirectly and thus D (£ m +i) only acts by substitution. 

The formula for R^ (|4.6|) follows after computing the resulting actions and rearranging 
the sums thanks to the re-summation formula provided by the contour integral: 



/ 7T 



Note that the parameters x a appearing in the contour integral (|4.12p are generic. □ 

Using the boundary-bulk decomposition of Proposition 13. 2[ one can now deduce 
from Proposition 14.11 the action of Q m (k) on arbitrary boundary states. 

Corollary 4.1 The action of Q m (k) on an arbitrary boundary state \ )& reads: 

m 

Qm («) I {A}f > & = £ £ ^ (Vxft) I {^} 7+ U {\} a _ ) b . (4.13) 
n=0 V x ; Vt- 
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The sum over partitions is defined as in Theorem \4-l\ and the coefficient can be 
expressed as 

K(Px,n) = E KA-Px,Vz) ^({£} 7+ ,{A CT } Q+ ), (4.14) 

o"i=± 

where S% ({^j™ , ) is the bulk function defined in ([4.8p . while lZ a {P\,V^) is the 
boundary dressing of (|4.7h : 



n \<k) n [/(Ab, ^) /(-A6, as; 
n {«(&) n w.eo n [Hh^m-h,^)) 



gw^wj (415) 



Here H? a y ({A} a+ ) and H B ({t;}~ /+ ) stand for the boundary-bulk coefficients f|3.8[) 
associated respectively to {\} a+ , {cr} a+ , and to {£} 7+ , { <J }7+ = {!)••• > !}■ 

Proof — The proof is a straightforward consequence of the boundary-bulk decompo- 
sition (|3.7p applied to Proposition 14. 1[ More precisely, expressing the boundary state 
|{A}^ )h in terms of the bulk states ({A""}^ ), and using (|4.3|) . we get 

m 

Qm («) i {A}f > 6 = E E E H w}(Wx) K (TV , I {ei 7+ u {n a _ ). 

n=0-p A ;-p f CTi=± 
l<i<iV 

We now use the fact that 

n /(-Af.A?) Hf ({A} , 

*W> ^ n TfTp^y -gfef «W <«>. - (^-), 

where ({£}7+ U {A}«„) is the boundary-bulk coefficient of |{^} 7+ U {A} Q _ )& 

in terms of U {A CT } a _ ). Note that the first factor of this product combines with 

the products over b £ a_ in the expression (|4.7j) of R(Vx",V^), and that the resulting 
factor, 

n im,K)f(-^K)] 

6 l a A n [f(K^a)f(-K^a)] ' 

ag7+ 

is actually independant of the value of o~i for i S a_. It enables us to reconstruct the 
boundary state |{£} 7+ U {A} Q _ )(,, with a coefficient which reduces to (I4.14p . □ 
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4.2 Action of local spin operators 



We list here the action of the local spin operators cr~, <t+ and E£f on bulk and 
boundary states. We omit the proofs since, although a little more technical, they 
parallel the one concerning the action of Q m (k). 

Proposition 4.2 The action of cr~ ; and <7+ on an arbitrary bulk state |{A}^) 
can be expressed as 

m— 1 

^ m |{A}f}=E E R-(V^,V^) |{£} 7+ U{A} Q _), 

m-1 N 

«iwf) = EE E ^f,n)im 7+ u{A} Q _), 

n=0 ci=l-p22 ^ 

m-1 N N+l 

a+|{A}f)= lim EE E E R n(n^) IUI 7+ U{A}f \{A} 5+ ), 

n=0 ci=l c 2 =l -p+ -p, 

in which the sums run over the following partitions 



{£}? = U} 7+ U {0 7 _ , I 7+ |=n + l, (4-16) 

Wf = {A} Q+ U {\} a _ , mtf* |«+|=n, (4.17) 



with 


1 7+ 


= n + 1 


with 


1 


1= n, 


with 


a + 


1= n, 


with 


a + 


|= n. 



: {Afc}i<fc<iv = {A} Q+ U {A} a , with | q+ |= n, (4.18) 

Vt ■ {h}i<k<N = {A} Q+ U {\} a , with I a+ |= n. (4.19) 

PFe afeo define the following partitions, associated respectively to (|4.18p and to (|4.19j) . 

= {A}f = {A} 5+ U {A} a _ , wift 5+ = a + U {<*}, (4.20) 

: {A}f +1 = {A} 5+ U {A} 5 _ , with 5+ = a+ U {<*, c 2 }. (4.21) 

r/ie coefficients R-(V^,V(.), R^(Vf,V^) and R+(V^,V^) are given as 

n smh(c a -e) 

^.P{)=^.^) gSt °n + smh(A a -e) ^^ A W'^ + ^> )' ( 4 - 22 ) 

Rf(Vf,V i )=R{Vf,V i ) sinhr? J] /(A a ,A Cl ) 

a6a+ 

n smh({„-{) 

AW--;) § -«AW,«K;l,{Ac,», (4.23) 
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Rt{Vl , n) = R{V+ , P e ) /(A C25 A C1 ) [] | sinhn [] /(A a , A C J | 



J] sinh(AAT + i -£ a + r)) 



aG7 + 



17 sinh(AAT + i - A a + rj) 



S n ( {£} 7+ , { A} a • Ajv+i , { A C1 , A C2 }) . (4.24) 



Here R(V\,V^) is given by (|4.7j) . and i/ie structure of the factor £n({£}i + , {A}™; £, {/u}^) 
is similar to (|4.8p . - 

II II sinh(£ 6 - A a + 7/) 

a=l 6=1 



^({e}? +1 ,{A}?;e,M?: 



n sinh(£ a - £ 6 ) J] sinh(A fe - A a 



x det n+ i 

where the matrix elements of M are obtained as 



M({A}?,U}" +1 ;C,M?) , (4.25) 



M{\}?, W 0*K)] iJfc = [^(A.M) ( Mi > (Oi )] iJfc . ( 4 - 2 6) 



with 



K(\j,{n}) = (1 - ^.n+ljll 



J /( A i>Mi)' 

in which 6{j denotes the Kronecker symbol and M K is defined as in ()4.9 



Using again the boundary-bulk decomposition, we are now in position to list the 
action of local spin operators on boundary states. 

Corollary 4.2 With the same notations as in Proposition ^-^ the action of o~~ n , 
and a^ t on an arbitrary boundary state |{A}^ )& takes the form 

m—l 

^|{A}f>6 = E E KnCPx'Te) l{0 7+ u{A} a _) 6 , 

m—l N 

^ 2 i{A}f), = EE E ^ 2 (^ 22 ^)im7 + u{A} Q _> b , 

n=0 ci=l-p22 ^ 

m-1 AT N+l 

^|{A}f>6 = , lim EEE E ^a + ^) IUI 7+ U{A}f \{A} 5+ ) 6 . 



n=0 ci=l C2=l -p+ p. 



The boundary coefficients 1Z , 1Z 22 and 1Z + have a structure similar to their corre- 
sponding bulk counterparts: 

n smh(c a -e) 

ae- l+ 



(Ti=± 



e-€m EI sinh(A a -C) 



Sn({A} a+ ,{£} 7+ ;£,0), 
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K 2 (r 2 x 2 ,V i )= ^{Vf^) smh ?? 11 f{X a , A c 



(Tj=± aga_| 



X , liln tt + „:,um FT S n(Wa + ,{£,}-y + ;t{X Cl }), 



n sinhfe-o 

«"Tm I] sinh(A a -£) 
K^'Pe) = E M^A>n) /(Ac 2 ,A^) njsinhr? J] /(A^A*) 

(Ti = ± 1=1 ^ 06a-|- 

EI [/(-^AT+1, £a) sinh(AAT+i - £ a + 7?)] 



X 



067+ 



-Sn({<!;} 7 ,{A CT } a ; Ajv+i,{A£}), 



II [/(-A 7V+1 ,A-) S mh(A7V +1 -A- + r ? )l "" VlS " " 



where lZ a is defined as in (|4.15|) and 5 n is £/ie bulk quantity (|4.25p . 



5 Correlation functions in the half-infinite chain 

We apply the results of the previous section to derive the expectation values of 
the generating function (Q m («)) of (<x^), and of (cr^ <7 m+1 ) in the ground state of the 
half-infinite chain. These are the boundary analogues of the results published in [49]. 



5.1 The generating function (Q m (ac)) 

Proposition 5.1 The generating function (Q m (k)) is obtained, in the thermody- 
namic limit M — ► +00, as the homogeneous limit of the quantity 

r+({?ir) Cd 

x det n [M K ({A},{z})] det n [V(\j,z k )], (5.1) 

in which M R is given by (14.91) . and W_ is the boundary dressing, 



n 2 



fT sinh ( Zj + £_ - r//2) 

W-({A}r,{4i 2 ) = ^ 

II sinh(A i+ £--r?/2) 

ni »i2 

n n sinh (^ + A a - v) 



U2 ni _ \ I J 

n sinh (z o6 - 77) sinh (A ab - 77) 
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of the bulk quantity 

m n 2 



W( {\}? , {z)T ) = = • (5.3) 



I! [ sinn ( z b - A a - 7?) sinh (z b - X a + 77) 
' =1 

n sinh (X ab - r/) n sinh (z ab + rf) 

a, 6=1 a, 6=1 

The contour of integration Co depends on the boundary magnetic field 

c f]-A;A[ur + (A) if o<f_<CA 

[ J —A ; A [ otherwise, 

where T± (z) stands for a small loop of index ±1 with respect to z. We recall that 
A = +00 for -1 < A < 1 and A = — mt/2 /or A > 1. 

Proof — Corollary 14,11 yields the action of Q m (k) on a boundary state. It is conve- 
nient to note that the coefficient TZa('Px,'Ps) flUISD can be rewritten as 

n 

M*^) = ( II °0 (^hr?)l^l J] W_({A CT } a+ ,{e} 7+ ) 
n sinh(£ a - £ 6 ) n sinh(A b - A a ) 

x ~ n n ^mI-k + „ } s.({a}„ + , {{ } 7+ ;(a}«_)-', (5.5) 

a£a + 6^7+ 

in which S<t({A} q+ , {£} 7+ ; {A} a _ ) is the function defined in (|3.16p . Then, using the 
reduced scalar product formula f)3. 15|) and absorbing the sums over partitions "Pg into 
auxiliary z integrals^], we obtain the former representation. 

Note that the contour contains T + (A) for large positive boundary field since we 
have to absorb the contribution coming from the replacement of the complex root A 
as explained in [4]. □ 



We refer the reader to [49] for technical details. 
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5.2 The ground state expectation value {(y^a m+l ) 

Using the same method as for the generating function (Q m ( K ))> we can also com- 
pute the ground state expectation value {^t a m+i)- It gives 



m—l 



sinh(^i + £_ - r//2) 



n=0 



n!(n + 1)! 



n+1 

m+1 EI f( z b,(,a) n 



n+1 j „ n+1 . 



n sinh(A b - £1) ^rj sinh(z fe - £ m+1 ) 



n+2 



6=1 



Q ft / ( A 6> £ a) smh(A 6 - £m+i) ^ sinh(z 6 - 6) 
6=1 

n+1 n+2 

[7 sinh(A fe - A„+i + rj) \[ sinh(A n+2 - A 6 + rj) 

^ - >v_ ( wr 2 , {zyr 1 ; 

n [sinh(z 6 - A n+ i + rj) sinh(A n+2 - z b + r/)] 

6=1 

n+2 

sinh(^i + A& - rj) 

6=1 



n+1 

n sinh(^i + z b - rj) 

6=1 



det n 



+i 



M({AK ! { 2 }r 1 ;Cm+l) 



x det n+2 Ci), ^(Aj, zi), . . . , y(\j,z n+1 )} . (5.6) 

In this expression, W_ denotes the boundary quantity (|5.2p . M( {A}™ , {z}™ +1 ; £ m +l) 
is a simplified notation for the matrix ik/^{A}-^ , {z}^^ >Cm+i> 

0) defined in (fl~26l) . 

Co is the contour (|5.4p . and denotes the following contour (^4-type contour): 

Ca= f ]-A + r/;A + 77[ur_(A), if - C/2 < £- < 0, 
\ ] —A + 77 ; A + 77 [ , otherwise. 

In the homogeneous limit, this results simplifies into 



(5.7) 



m— 1 • 1 > ,> n+1 , „ n+1 „ 

n_0 r+(n/2) fe - 1 C D k ~ l C A 



n+1 

n 

a=l 



sinh(2; a + rj/2) 
sinh(z a - 77/2) 



n 

0=1 



n+2 

J! sinh(A 6 - t?/2) 
6=1 



sinh(A a - 77/2) 
sinh(A a + r ? /2)J -+ 1 ■ , 
6=1 



n+1 n+2 

]1 sinh(A fe - A n+ i + rj) \[ sinh(A n+2 - A& + 77) 
x ^ ^ W_ ( {A}" +2 , {z}^ + 1 ) 



n+1 



fj [sinh(z b - A n+ i + rj) sinh(A„ +2 - z b + rj) 
6=1 



x det 



n+1 



M({A}™ , {z} n + l ; 77/2) det n+2 [*(A ia 77/2), * (A,-, Zi), . . . , *(A i; z n+1 )] 
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6 An alternative resummation 

6.1 Bulk type resumations 

We have obtained in the previous sections a series representation for the generating 
function. It happens, just as in the bulk case [52], that it is also possible to derive a 
totally different representation for (Q m (k)). The latter is based on a re-summation 
of its expansion with respect to elementary blocks : 



(Q m («)> = (II 1 + ^ 22 )> = E k ' f - . 



1=1 



s=0 



where 



1 



s\ (to — s)l 



2, i = l...s. 



s + l...m, 



(6.1) 



(6.2) 



and S m is the group of permutations of m elements. These elementary blocks were 
computed in [4]. They can be written as multiple integrals in the half-infinite size 
limit: 

^r' 1 • ■ -e^-) = (-ir- s ffidxJ n dA, - ^MhM — 

c D i=1 L i=s+1 11 [sinh^- sinh (^ - rj) } 

i<j 



]1 sinh (Xi + £j - rf) 

x 7= 7 

n sinh (Xij - rj) sinh (Ay - rj) 

i>j 

TT sinh ^ + ~ r] / 2 ^ TT 

J-J-sinh(Ai + C--r?/2) 

1=1 v % s " ; p=s + l 

The indices i p are defined by 



n 

P =i 



i„ — 1 



| I sinh (£,■ - A p ) J J sinh (£,■ -X p -rf) 

j=l j=i P +l 



H sinh - A p ) H sinh - X p + rj) 

j=l j=i P +l 



(6.3) 



{i : 1 < i < m, e- = 2} = {n < ■ ■ ■ < i s } , 
{i : 1 < z < m, ej = 1} = {i s+ i > • • • > i m } . 



(6.4) 



For simplicity, we consider from now on the massless regime (although all what 
follows can be performed in the massive regime as well). In that case, rf = —i( and 
the contours of integration Cd and Ca depend on the boundary magnetic field h- as 
follows: 



range of 


D — contour 


A — contour 


C/2<|e-|<vr/2 


= R 


C A = K - K 


C/2 > £_ > 


c D = MUr+ (A) 


C A = K - «C 


-C/2 < £- < 


C D = R 


C A = {M-iC}Ur- (A) 



(6.5) 



We recall that A = r//2 — and that r±(z) is a small loop around z of index ±1. 
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We now perform a change of variables in the A-type contours: = — iQ. 
Moreover, we shift the inhomogeneities around zero Si = £j + i£/2, and define 



cii = 3/2 - a 



1/2 (q = 1) for A-type, 
-1/2 (ei = 2) forD-type. 



(6.6) 



This gives 



..El 



r(m) ^7r(m) \ 



(_1)W J d s\ J d ™s\ 
Cd C a 



det r 



n*(Mi; 



n 



sinh (5fc - A^y) + ia,r(j)C) sinh [Sj - A 



7r(fc) — ia ir(k) 



)0 



sinh (A 7r 0) 7r ( fe ) - mj^x^C) sinh (A 7r ( J -) 7r ( fc ) - ia 7r (j) 7r (fc)C) 



Here 



*(A,5)=*(A,5-iC/2) 



p(A-j^-p(A+j) 
2 sinh 25 



-A;A[ur_ (A-i)) if-C/2<£_<0, 
-A ; A [ otherwise. 



(6.7) 

(6.8) 
(6.9) 



and (— l)' 71 "' is the signature of the permutation. 

One can compute the sum over permutations (|6.2[) just as in the bulk case [52]. 
It leads to the following integral representation for F s : 

Proposition 6.1 (Bulk- type resummation) The generating function of the spin 
correlation function (Q m (k)) can be expressed as 



(6.10) 



with 



s\ [m — s)\ 



det T 



d s A / d m ~ s A 



s=0 



n 



sinh (£_ + Sj ) 



Cd C A 



J] S (6i,Sj) A = -L sinh (Xj + £_ - idj() 



Kj 



x 6 S ({A}) Z m ({A} , {5}) J] sinh (A,- + <S fc - iC«i) • (6.11) 

j,k=i 

Here Z m ({A} , {£}) stands for the partition function of the six-vertex model with do- 
main wall boundary conditions: 



ns(A i -<5 fe ,zC/2) 
m } 5 { }) ]^[ sinh Aj/c sinh <5fcj m 



1 



s(Aj - S h ,i(/2) 



(6.12) 



20 



while 



».({*» =n n - ; inh ui - rt 



(6.13) 



Similar representations can be obtained for other correlation functions. Here we 
give only two important examples: the local density of energy and the (^m+i a i) 
two-point function. 

The local density of energy 



Em = {<T m <T X m+ l + « +1 + ~ 1)> 

can be written as a sum of m terms 

m— 1 



s=0 



E m — ^2 Esi 
each of them containing m + 1 integrals 

_\ _ g) , 1 d s A J d ro — X A I d\ m J dA m+1 ^^ 



(6.14) 



(6.15) 



det r 



*(Ai,^-) 



C D Ca Cd Ca ' l< 'j 

x#s ({Ai, . . . , A m _i}) Z m ({Ai, . . . , A m _i} , {<5i, . . . , <5 m _i}) 

nsinh (£_ + 5j) '^y 1 sinh(Aj — 5k — ia-jC) sinh(Afc — 5j + ia^C) 
^h7I7T7 -in.;rt 11 11 



3 

m+1 



L sinh(A j + £_ - iajO 



j=m k=l 



sinh(AjA: + iajkC) sinh(Ajfc - ia jk Q 



n. , . c . ).\ ^(Am, A m+ i, <5 m , <5 m+ i) 
smh (A,- + Ofc — m,<) . . : — v . , ; r, 
h=1 { J 3 ' smh(A m+1 - A m ) smh(A m + A m+1 ) 



(6.16) 



where 



(p(\m,\m+i,5 m ,5 m+1 ) =sinh(A m - 5 m+1 + i^) sinh(A m+ i - 5 m+1 - S) 

c C 

+ sinh(A m -5 m - i-) sinh(A m+ i -5 m + i-) 

- cos C sinh(A m - 5 m+1 - i^) sinh(A m+ i - 5 m - S) 

-cosCsinh(A m - 5 m + sinh(A m+ i - 5 m+1 +^)- 

A similar representation can be obtained for the two-point function (o"^ +1 cjf ), 
namely: 

m— 1 
s=0 
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and every term contains m + 1 integrals 



G, 



(-1) 



m—l 



s\ (m — 1 — s)! 



det 



S A- 



m+1 



d s+ ! A 

Cd c a 

Xds ({^2, • • • , A m }) Z m ({A2, • • • , A m } , {62, ■ ■ ■ , S m }) 



n * c^ii 

i<j 



m+1 



sinh (£_ + (5j) sinh(Aj — 5\ — iaj() sinh(Aj 
sinh (Aj + £_ — idj() 



$m+l +icij() 



m+1 m 

x 11 smhiXj + Sk- iajC) J] II 1M - , • 1 ,t - , 

J= r; m +i k=2 smh (V + lOjfcC) smh(A ifc - ia jk Q 



sinh(Aj — 5k — iajC) 



x ■ 



sinh(Ai — 61 — ii) sinh(A m+ i — Si + ii) 



sinh(A m+ i - Ai) sinh(Ai + X m+ i) 
6.2 Boundary type resumations 



(6.18) 



It is important to note that the function Z m appearing in the representations 
()6.1ip . (|6,16p and (|6,18p is the bulk partition function represented in terms of the 
Izergin determinant [71]. In the boundary case one can symmetrize the integrand 
even further by writing it in a form invariant under the reversal of the parameters 
A and finally rewrite the result in terms of the boundary partition function and the 
Tsuchiya determinant. The integration contours in (16. lip are not invariant under 
the transformation A —* — A. We thus deform the contours until we obtain a reversal 
invariant contour. As we do not cross any pole of the integrand, the result remains 
unchanged. Actually we can even pick the contours so as to integrate only over 
one contour C, although this is not necessary. This contour C is defined as follows 
according to the value of the boundary field h-\ 



C/2<|£_ |<vr/2 


C = R 


c/2 > i- > 


c = r u r+(A)U r_(-A) 


-C/2 < £_ < 


c = r u r + (ic + A)U r_(-<-A) 



(6.19) 



We extract the totally even part of the integrand appearing in (16. 7|) according to 



J dxf(x) = ±J2 j dx f( x<T ) . 



x = a x . 



(6.20) 



We get 



s!(m - s)!2 r ' 



/ 



det r 



d m A 



*(Ai,*,-)J sinh (£_+<*. 



Kj 



?= L sinh(A i ,^_ - iajC) 
x6 s ({A}) H s ({\},{5}) 



(6.21) 
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Here 

Q s ({A}) = TT _*( A i' A fe) _ ( (6 . 22) 

j>k 5 ( A J' fc ' m jkO 5 {^jk,iaj k Q 

and the sums over negations have been absorbed into H s ({A} , {5}): 

m sinh ( A!L + iajkC ) 

h s ({A} , {5}) = e n & sinh ( a j - - ] n — — - 

o-i=±j=i j>k sinh ( A jfc 

m 

x 11 sinh (AJ + <J fc -ia J -C) Z m ({A CT } , {£}) . (6.23) 

Equation (16.231) implies that H s ({A} , {8}) is a symmetric function of the parameters 
A and of the parameters 5. Moreover, Q 2 ( m - l )^j fj s ({A} , {8}) is a polynomial in each 
of the e 2Aj variables of degree 2(m — 1). Finally, it is a matter of straightforward 
computations to check that H s satisfies the reduction properties: 

H s Ui=±(5i-ioiC) (i^iLi ; {Sk}T=i) = ±H S ({AJ™ 2 ; {4}£L 2 ) 

m 

x sinh(2(5i - ia-yQ) sinh {8\ - £_) TJs (Aj, <5i +ioiC)s(5i -2zaiC,4)- (6.24) 

3=2 

These are the reduction properties of Z m ({A} , {5}), the partition function of the six- 
vertex model with reflecting ends [72] . Supplementing this result with the equality of 
the two functions at m = 1, we obtain that H s ({A} ; {8}) is s-independent and equal 
to Z m ({A} , {8}). Hence, we have the following result: 

Proposition 6.2 (Boundary-type resummation) The generating function of the 
spin correlation function {Q m (n)} can be expressed as 

m 
s=0 

with 



F s 



1 f det m |tt (Aj, 8j)j ™ S jnh(^ + ^ 

-s)\2™J U 8 ^ 8^ ij- sinh(A J ,e--« 



s!(m-s)!2 m 7 Il-s(<Mi) M sinh(A i ,^_ - iajQ 

C i<j 3 

x 6 S ({A}) Z m ({\},{8}) . (6.26) 

Here 

O s ({A}) = TT -1%^ _ , (6.27) 
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and 



Zm({\},{6}) 



II [s(A j ,5 fc + </2)5(A J ,5 fc -</2)" 
j,k=i 



J] [s(\i,\j) s{5j,5i 



x JT [sinh2Aj sinh((5j — £-)] det r 
i=i 



1 



s(A J ,4 + iC/2) 5 (A J , ( 5 fc -iC/2) 



(6.28) 



Let us finally mention that the sum over s in (|6.2[) can be rewritten as a single 
integral over auxiliary z variables: 



(Qm («)} 



1 



m!2 r 



d m A 



d m z 

2 m ({X},{5}) det m ^(Ai,^) 



U r+(«$) 



(2ivr) r 



n [5(A j ,A fc )s(A jfc ,«C)s(A jfc ,iC) 

j>k 



X H s(6i,5 k ) 5 (A p ,e_ + 



fc>Z 



where 



sinh 2z K -*(*+</2)/£ 



(6.29) 



(6.30) 



a(MC/2) 

This boundary-type resummation yields an integrand not only symmetric in {A} but 
also invariant under a reversal of any integration variable A. These properties allow to 
compute completely the so called emptiness formation probability at A = 1/2 and for 
vanishing boundary magnetic fields. It also allows to obtain the leading asymptotics 
of this quantity at the free fermion point [75]. 



7 The free fermion point 

7.1 Local magnetization at distance m 

The first (and the most important) application of the re-summation methods given 
above is the magnetization profile. This one-point function 

(a z J = l-2D m d K (Q m («)) | K= i, 

can be computed at the free fermion point by using the two different types of re- 
summations for the generation function. We give here both derivations. 

7.1.1 First method 

In the free fermion point, the n th term of the series ()5.ip behaves as (k — l) n . 
Thus, after taking the n derivative and sending k to 1, only the n = 1 term survives. 
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At C = vr/2, we have 
<Gm(«)> = 



Thus©, for m > 2, 



■11 



\ {2nf n (n!) 2 7 7 lllltanh^-e.) 

n=0 v ' v ' n r , ftl , a=l 6=1 



nf sinh (zj + £_ + nr/4) 
^ x \sinh(A J - + £_ +Z7T/4) 



sinh 2A, 



xdet r 



sinh (Aj — Zfc) 





1 


det n 


_s(Xj,z k )_ 





dX 



sinh (A - - 27r/4) [tanh (A + in/A)] 



2(m-l) 



,t j sinh (A + £- + W 4 ) cosh 2 (A + ivr/4) 



Computing, if it exists, the residue at A we get, for m > 2, 



(7.1) 



(7.2) 



«> = -28 (/»_-!) 



+ 



(-1)' 



7T 



dA 



/i_ + z tanh (A + i7r/4) [tanh (A + nr/4)] 



2(m-l) 



(7.3) 



1 + ih- tanh (A + in/4) cosh 2 (A + in/ 4) 

K 

where (x) is the Heaviside step function. The standard A = change of variables, 

e ip = -tanh(A- wr/4) , (7.4) 

yields 

7T 

hi - 1 



-2Q(h- 



1) 



h 2 * 



-JT f dp e -2i(m-i) P e ; p + ih. 
n J e tp — ih^ 



(7.5) 



Thus (cr^j) displays Friedel type oscillations induced by the boundary. Moreover it 
decays as 1/m when m — > +oo: 



2(-l) m h. 



nm 



hl + l 



+ -5 , 



tu- 



rn » 1. 



(7.6) 



Here we recover the results of [69], since we have /i_ = \/2a_ in Bilstein's nota- 
tions. When | h- \— > oo we conclude from (|7.5|) that the first site is totally decoupled 
from the others as {(T^ n ) m >2 goes to its bulk average value 0. Actually in this limit 
the model is in correspondence with a Kondo model with a spin 1/2 impurity [29]. 



We do not give (of) as it corresponds to d K (Q\ («)) without taking the lattice derivative. 
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But in this case the impurity is completely screened, and the overall magnetization 
in zero. 

We also recover from (17. 5h . just as expected from the spin reversal symmetry, that 
(°m) = when the boundary field vanishes. Actually this observation holds for all A 
as inferred from the structure of the monodromy matrix (12. 6h on the first site. When 
one sets £i = rj/2 and = then it acts as a diagonal matrix on the first site, a sign 
of the claimed decoupling. 



7.1.2 Second method 



Starting from the re-summation formula (|6.1ip of Proposition 16.11 we implement 
the simplification due to £ = it/ 2. If we perform the change of variables 



e ip = -tanh(A-ivr/4) 
in (|6.1ip at £ = tt/2 then we arrive at 



(7.7) 



m(m — 1) „ 

s\ (m — s)\ 1A - 

s m 



2vr 



e iPj _|_ e -Wj 



n 



2ir 



e «Pj _|_ e ~Wj 



C A 



X U II [(e- ipfe -e^)(sm Pi + sin^)] [] [ 

k=l j=s+l j,k=l 

j>k 



) (ship, 



smpk) 



(7.8) 



x Yl [(e^ -e^)(sinp fc -sin Pi )]. 

j,k=s+l 

j>k 

The contours of integration are 

Ca = ]0;tt[, h->-l, C A = ] ; vr [ U T_ (e" ip = -ihJ) , /i_ < -1, (7.9) 
Ca = ]-tt;0[, /i_ > -1, C A = ] -vr;0[ur + (e ip = -ihJ) , /i_ < -1, (7.10) 
C D = ]0;n[, h-<l, C D = ] ; vr [ U T + (e ip = -ih-) , h- > 1. (7.11) 



Once we introduce the function 

<?«(p) 



k peC D , 
l psCa, 



(7.12) 



we can re-sum the terms F s into a single m-fold integral for {Q m (k)}: 

m(m—l) 

(Qm («)) 



(2i)- 



mi 



CaUCd 
m 



n 

i=i 



2;r Jk( ^ (.'/-, - ;/,_ 



n 

j>fc>i 



) ( sin ^ 



smpfc) 



(7.13) 
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One can then express the generating function as a single determinant 

(Qm («)) = det m p7(K)], (7.14) 

If Jpk _ (_l)k -ipk 

U jk (n) = - J d P 9 K (p)e-^-V J_> h _ . (7.15) 

C A <JC D 

To simplify this result we add to each row of U(k) the next one multiplied by i/t_: 
(Q m ( K ))=det m U^(K), 

U^\k) = i- y dp fl K (p) (e ip ( fc -^ - (-l)*e-^ + *)) , j<m, 

C A UC D 

It is easy to see that <3m(l) = 1- Computing the first derivative of the generating 
function one recovers the result already obtained from the series (IT. 5j) : 

/\ -Mprm-iie-* + ifc- 



(-l) m /■ 



e^ - ih. 

c D 

h - ~ ^ , (-!) m L„ „-2ip(m-i) e ~'" - 



-2^B - 1) + ^ / dp (7.16) 



o 



7.2 Local density of energy 

The local density of energy is another interesting quantity [61] that one can eval- 
uate for the XXO chain: 

E m = (cr m cr m+1 + C7 y m a v m+1 ) . (7.17) 

Starting from (|6.16p and using the same technique as in the previous sub-section one 
easily obtains the following representation for the density of energy: 

E m = -2ij^- y g (e^+s) + l) det m+1 [V(p,q)] . (7.18) 
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The entries of V (p, q) read 
Vjkip, 0) = ^ J dp' (e^( fc -J) - (-l) fc e -^'0+ fc )) = 5 jk , j<m-l, 



C A UC D 

1 /• Jv'k ( -\\k -ip'k /u \k-m+l 

V m - lk ( P , q ) = ±- Uj/e-^-Q * -J- 1 ^ 

c A uc D 



— ) 0(ife-m+l), 



3 ipfe _ (_i)k e -»i* 



e^ - 



Wi fc (P,9) = e- J,?m - 



e*<? - ih- 

Finally, the integrals over Ca can be represented as 

/"/-/• 

Ca Ca^JCjj Co 
Accordingly, E m reduces to a sum of two 3x3 determinants: 



(7.19) 



(7.20) 



E m = -2i 



where 



1 ih- -h 2 _ 

F(m — l,m—l) F(m — l,m) F(m—l,m+l 
1 ih- 

1 ih- -hi 

—2i F(m,m—1) F(m,m) F(m,m + 1 
1 



If . . Jpk _ (-l)k e - 

'0\*) = ^ Jdpe-™- 

C D 



e'P - ih. 



(7.21) 



(7.22) 



The computation of these determinants yields 

E m = -2i (F(m, m) - F(m - 1, m + 1)) - 2h- (F(m, m - 1) - F(m - 1, m)) , 

(7.23) 

or, more explicitly, 

Em = -- + - (-l) m / dp e-^ 2 ™" 1 ) e ~J + f- . (7.24) 
7r 7r y e l P — ih- 

C D 

The constant term reproduces the bulk result. The influence of the boundary 
appears in the oscillating term. In the m — ► oo limit the local density of energy 
behaves as 



4 2 1 — h 2 

E m = -- + —— (-l) m 7— ry- 
7r 7rm 1 + h- 



10 A 



(7.25) 
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7.2.1 Two-point function (o"^, 1 a 1 ) 

The preceding method can be successfully applied to compute other types of two- 
point functions (like boundary-bulk two point functions), here we give the example 

In the free fermion point, after the usual change of variables and some straightfor- 
ward but tedious calculations, one obtains from (16.181) a simple determinant formula 
for this object. 

(^m+i^D = -idet m+1 V+-, (7.26) 

(7T 2tt\ 
J - J\ dp ( e *K*-i-D - (_i)* e HKj+i+*)) , j < m — 1, 

1 [■ ipk _ I _-\\k „-ipk 



2vr J e l P - ih- 

1 f pipk _ (_-\\k p -ipk 

Computing the integrals in the first m — 1 rows and using the fact that sum of the 
last two rows is 5k ; m+i we reduce this representation to a determinant of a m x m 
matrix 

("™ + l<'r)= J 2^-det m F + -, (7.28) 



(j + l) 2 - k 2 



ipfc 



C>-=jdpe-^ e ,;_,;_ , (7.29) 

Co 

This determinant can be computed for any value of m. However the result is quite 
different for m odd or even. The details of the computation are given in Appendix 
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B, here we give only the final result for the two point function 



)2a-l 



2a-l 



n 



r(j) \ r(a-i)r 3 (a + i) 



\ r(j + \) J r(o)r(a + i) 

-jp(2a-l) 

ih- 



f ( \ e _ip(2 

C J dpP 2a _! ^COSp) — - 



1^ T(j) \ Via - §)r»i« + & 



r(o)r(a + i) 



(7.30) 



dq cos q j dpP 2a (^cos(g - p)j - . p _ — , (7.31 



where P m (x) are Legendre polynomials. 

Asymptotic analysis of these expression yields the following leading behavior of 



KUi^n =(-ir^_)m-f 



2 ri r 00 ^ 

ex P It / j 



vr(l + /i 2 ) "" r I 4 



1 



cosh t 



(7.32) 
(7.33) 



8 Conclusion 



In this article we have obtained different types of physical correlation functions 
of the open XX Z chain from re-summations of the multiple integrals derived in [4] 
for the elementary blocks. At the free-fermion point, we were able to use these 
representations to derive explicit results such as the formula for the density of energy 
profiles, a quantity arising in the study of quantum entanglement in spin chains [61]. 

Just as in the bulk case, the question concerning the asymptotic behavior of the 
correlation functions outside of the free-fermion point naturally arises. The problem 
is of the same order of difficulty as in the bulk model. Indeed, the multiple integrals 
differ from their bulk counterparts only by factors due to the Z 2 symmetry A — > — A 
and the presence of boundary fields. 

One could also wonder if it would be possible to tell something about the dynam- 
ical or temperature correlation functions. It seems that this generalization is highly 
non-trivial. 

Finally, we would like to stress that our expressions also simplify at other partic- 
ular points such as A = 1/2. For instance, when A = 1/2, one can already compute 
completely the so-called emptiness formation probability when h- = [75]. 
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Appendices 

A Asymptotic of the two-point function (cr^ +1 crf ) 

In the last section we obtained a determinant representation (|7.29|) for the two- 
point function This determinant can be computed for any value of m. 
However the results are quite different for m odd or even. 

If m is odd: m = 2a — 1, the result can be written in the following form 



det m y H 



TT 



n 



i r(i + h: 



r(q-i)r 3 (a + i) 

T(a)T(a + l) 



E 



r(a-6+i)r(a + 6-i) 



r(a-6+l)r(a + 6) 

ip(26-l) + e -ip(2b-l) 



b=l 



- ih- 



(A.l) 



If m is even: m = 2a, the result is quite similar but there is a very important 
difference: 



det m y+" =2 m+ V m - 3 



n 



r(a + |)r 3 (a + i; 



i r (j + i) J r(a)r(a + i) 

T(a-6+i)r(a + 6+i) 
J (6 + i)(6 - i)"r(o-6+l)r(a + 6+l) 



*E 

6=1 

TT 

x J dpe 



-2ipa 



2ipb 



-2ipb 



e l P - ih- 



(A.2) 



Asymptotic analysis of the prefactors in (jA.ip and (|A.2|) is rather simple, namely: 



'2a-l 



n 



r(i) 



r ( a-i)r 3 (a + i; 
r(o)r(a + i) 



V^TT _1 

: m 4 exp 



2a 

n 



r(a)r(o + i) 



^/2vr 3 

2'm+l 





r°° di 




io T 


+ 2- 


) 


1) 






r*f 


{I 


'o * . 



-4t 



cosh 2 t 



-to 



1 



cosh 2 t 



i + o(- 



m 



i + oH 



(A.3) 



(A.4) 
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For m odd the sum in (lA.ljl can be rewritten as follows: 

" T{a-b+\)T{a + b- \ 
^ T(a-b+l)T(a + b) 



2 ' I e -2ip(a-b) , e -2ip(a+b-l) 



«=0 



r(/ + i)r(2a-z) 



(A.5) 



and can be represented in terms of the Legendre polynomials P m (cosp) 

2a-l n/, , hn/n. t In t-,/o hrtl 



rq + i)r(2a-/-i) c _ 2 ^ 



Z=0 



r(z + i)r(2o-z) 



r(2 °-' )r ® 2Fl (i,l-2„;|-2«;e-*) 



r(2a) '~ iV 2' 
=vre- i f m P m (cosp) 



Using Laplace asymptotic formula, 



P m (cosp) 



\7rm sinp 



2 \i 

I cos 



p I m + 



1 



7T 



+ o 



3 

ma 



for the remaining integral we obtain the following leading term 

p —ipm 



dpP m (cosp) 



e l P - ih. 



o 



TT 



m(l + /r ) 



l + 0( 



(A.6) 



(A.7) 



(A.8) 



Assembling all the contributions we obtain the following leading term for the two- 
point function (for m odd): 



■iy 



7T(l + h 2 _) 





r°°dt 




Jo t 



-It 



1 



cosh 2 t 



m~i (1 + 0{ 



1 



m 
(A.9) 

The same result holds for m even, but the derivation is a little bit more tricky. 
The sum in (|A.2|) can be once again rewritten in a more simple way 



T{a-b+\){h_)T(a + b + 



' x (b+\){b-\) T(a-b+l)T(a + b+l) 



dpe 



-2ipa 
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4e 



1 1 ^ r(a-b+ ±)r(a + 6 + ±) 



f>=— a 



6 + 1 6 _ i / r(o - 6 + l)r(a + 6+1) 



7T 
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dp 



2ipb g— 2ip6 

e*? - ih- 

e 2ip(b-a) 

e*P - ih- 



ys r(a-6+±)r(a + 6+±) 
r(a-6+l)r(a + 6+l) 



<e- 2iqb cosq / dp 



2ip(b—a) 

e*P - ih- 



■ J dq cos q J 



in \ dq cosq / dp P m (cos(q — p)) 



-imp 



e}P - ih. 



(A.10) 
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where we introduced an additional integral to be able to express the result once again 
in terms of the Legendre polynomials. Asymptotic analysis of these integrals gives 



3 



/I" e -imp / 7T \ f 2% ( 1 

dq cosq / dp P m {cos(q - p))-r- — = - ( — ) - l + 0(-=] 







(A.ll) 

and it leads once again to the same leading term (|A.9|) for the two-point function. 
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